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Abstract 

We study an extension of Axelrod's model for social influence, 
in which cultural drift is represented as random perturbations, while 
mass media are introduced by means of an external field. In this sce- 
nario, we investigate how the modular structure of social networks 
affects the propagation of mass media messages across the society. 
The community structure of social networks is represented by coupled 
random networks, in which two random graphs are connected by inter- 
community links. Considering inhomogeneous mass media fields, we 
study the conditions for successful message spreading and find a novel 
phase diagram in the multidimensional parameter space. These find- 
ings show that social modularity effects are of paramount importance 
in order to design successful, cost-effective advertising campaigns. 



1 Introduction 

Over the last few years, statistical physics has increasingly contributed use- 
ful tools and valuable insight into many emerging interdisciplinary fields of 
science pfl El E]- In particular, many efforts have focused recently on the 
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mathematical modeling of a rich variety of social phenomena, such as social 
influence and self-organization, cooperation, opinion formation and spread- 
ing, evolution of social structures, etc (see e.g. [5l [6^ [7| [8t [9j [TTTt [TT1 [T2t fTai 

in da Hum dEi ehj Eo]). 

In this context, an agent-based model for social influence, originally pro- 
posed by Axelrod [2H [22] to address the formation of cultural domains, 
has been extensively studied within the sociophysics community (see e.g. 
[231 EH ESI [261 EH [2HI EHl EU]). In Axelrod's model, culture is defined by 
the set of cultural attributes (such as language, art, technical standards, and 
social norms [22]) subject to social influence. The cultural state of an indi- 
vidual is given by their set of specific traits, which are capable of changing 
due to interactions with their acquaintances. In the original formulation, the 
individuals are located at the nodes of a regular lattice and the interactions 
are assumed to take place between lattice neighbors. Social influence is de- 
fined by a simple local dynamics, which is assumed to satisfy the following 
two properties: (a) social interaction more likely takes place between individ- 
uals that share some of their cultural traits; (b) as a result of the interaction, 
their cultural similarity is increased. 

Earlier investigations showed that the model undergoes a phase transition 
separating an ordered (culturally polarized) phase from a disordered (cultur- 
ally fragmented) one, which was found to depend on the number of different 
cultural traits available |23j. The critical behavior of the model was also stud- 
ied in different complex network topologies, such as small-world and scale- free 
networks [21]. These studies considered, however, zero-temperature dynam- 
ics that neglected the effect of fluctuations. Following Axelrod's original idea 
of incorporating random perturbations to describe the effect of cultural drift 
[21] , noise was later added to the dynamics of the system |25j. With the 
inclusion of this new ingredient, the disordered multicultural configurations 
were found to be metastable states that could eventually decay towards or- 
dered stable configurations, depending on the competition between the noise 
rate and the characteristic time for the relaxation of perturbations. 

Very recently, other extensions of the model were proposed, in which the 
role of mass media was investigated within different scenarios. Neglecting 
random fluctuations, some studies considered in detail the role of external 
[26] and autonomous local or global fields [27] . Another recent investigation 
focused on the interplay and competition between cultural drift and mass 
media effects [29] . Adopting a mass media coupling capable of affecting 
the cultural traits of any individual in the society (including those who do 
not share any features with the external message), it was shown that the 
external field can induce cultural ordering and reproduce the trend of actual 
advertising campaign data. 
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In a related context, recent investigations addressed the role played by 
the underlying topology of complex substrates on the dynamical and critical 
behavior of the models defined on them. The effects of some structural 
properties that characterize disordered substrates, such as the small-world 
effect, the degree distribution, the degree-degree correlations, and the local 
clustering, were extensively studied [311 E2J EHl Ell ES] • 

Furthermore, the property of community structure, or large-scale cluster- 
ing, appears to be common to many real-world networks and is nowadays 
being subject of intense research efforts. In many social networks, as e.g. the 
well known karate club study of Zachary [36], the United States House of 
Representatives [37], scientific co-authorships and mobile phone call records 
[3S] . well defined modular structures were observed. However, the effects of 
community structure on models of sociophysical interest have received so far 
little attention. Lambiotte and Ausloos [39] have very recently considered the 
effect of communities on the Majority Rule model by means of the so-called 
coupled random networks, a mixture of two random communities, in which a 
parameter v controls the degree of intercommunity links relative to that of 
intracommunity connections [3H] • Depending on v and on a noise parameter, 
a diagram with three distinct phases is obtained: a disordered phase, where 
no collective phenomena takes place; an ordered, symmetric phase, where 
both communities share the same average state; and an ordered, asymmetric 
phase, in which different communities reach different states. 

The aim of this work is to investigate effects arising from the character- 
istic modular structure of social networks in the propagation of mass media 
messages across the society. To this end, we focus on the extension of Ax- 
elrod's model proposed in Ref. [29], which includes effects of mass media 
and cultural drift, using coupled random networks for the substrate. In the 
absence of external messages, a phase diagram with three phases is found, 
qualitatively analogous to that observed in Ref. [39] for the majority rule 
model. Then, we assume that an inhomogeneous mass media field affects 
one of the communities and study the system's response to the spreading of 
the message. Incorporating the intensity of the mass media field as an ad- 
ditional parameter, several new phases are observed to emerge, thus leading 
to a very rich, novel phase diagram in the multidimensional space of model 
parameters. 

This paper is organized as follows: in Section 2, details on the model and 
the simulation method are given; Section 3 is devoted to the presentation 
and discussion of the results, while Section 4 contains the conclusions. 
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2 The model and the simulation method 



In order to represent the community structure observed in social networks, 
we consider a substrate topology consisting of two coupled random networks 
(CRN). These structures were first proposed in Ref. [36] to carry out compar- 
ative tests of different methods for community detection in complex networks. 
We assume that a system of TV nodes is divided into two communities (A and 
B) of equal size. A CRN configuration is built by adding intracommunity 
links between pairs of nodes that belong to the same community, as well as 
intercommunity links between pairs of nodes that belong to different commu- 
nities. Considering all possible node pairs, intracommunity links are added 
with probability pj nt , while intercommunity connections are added with prob- 
ability p ext . On average, a node is thus connected to k int = pi nt (N/2 — 1) 
neighbors inside the same community and to k ex t = p ex tN/2 nodes that be- 
long to a different community. For the sake of simplicity, we fix ki nt = 4 
and tune the intercommunity connectedness by means of a single parameter, 
namely v = k ext /ki nt ~ Pext/ Pint- Notice that the v — > 1 limit corresponds to 
a single random graph lacking modular features, while v <C 1 is the case in 
which well defined communities are sparsely connected with each other. 

The nodes of the system are labeled with an index % (1 < i < N) and 
represent individuals subject to interactions with their neighbors (i.e. other 
individuals directly linked to him/her by either intra- or inter-community 
connections), as well as with an externally broadcast mass media message. 
According to Axelrod's model, the cultural state of the i— th individual is 
described by the integer vector a- x = (<7ji,Ci2> •••jCif), where 1 < cr^/ < q. 
The dimension of this vector, F, defines the number of cultural attributes, 
while q corresponds to the number of different cultural traits per attribute. 
Initially, the specific traits for each individual are assigned randomly with a 
uniform distribution. Similarly, the mass media cultural message is modeled 
by a constant integer vector /z = (/ii, [12, Hf), which can be chosen as 
H = (1, 1, 1) without loss of generality. The intensity of the mass media 
message relative to the local interactions between neighboring individuals is 
controlled by the parameter M (0 < M < 1). Moreover, the parameter r 
(0 < r < 1) is introduced to represent the noise rate [24J. 

Since the main focus of this work is on mass media spreading phenom- 
ena under the influence of an underlying modular substrate, we will consider 
inhomogeneous mass media affecting only individuals that belong to the com- 
munity A. The model dynamics is defined by iterating a sequence of rules, 
as follows: (1) an individual is selected at random; (2a) if the individual be- 
longs to the community A, he/she interacts with the mass media field with 
probability M, while he/she interacts with a randomly chosen neighbor with 
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probability (1-M); (2b) if the individual belongs to the community B, he/she 
interacts with a randomly chosen neighbor; (3) with probability r, a random 
single-feature perturbation is performed. 

The interaction between the %— th and j— th individuals is governed by 
their cultural overlap, = ^2f = i8<r i f,<r j f/F, where 5ki is the Kronecker 
delta. With probability C^, the result of the interaction is that of increasing 
their similarity: one chooses at random one of the attributes on which they 
differ (i.e., such that o"i/ 7^ o"j/) and sets them equal by changing one of 
their traits. Naturally, if CV, = 1, the cultural states of both individuals are 
already identical, and the interaction leaves them unchanged. 

The interaction between the i— th individual and the mass media field is 
governed by the overlap term C{m = E/=i ^<t</,a*/ + -QA-^ 1 + !)• Analogously 
to the precedent case, Ciu is the probability that, as a result of the interac- 
tion, the individual changes one of the traits that differ from the message by 
setting it equal to the message's trait. Again, if du = 1, the cultural state 
of the individual is already identical to the mass media message, and the 
interaction leaves it unchanged. Notice that Ciu > 0; thus, the mass media 
coupling used here is capable of affecting the cultural traits of any individual 
within community A, including those who do not share any features with the 
external message. 

As regards the perturbations introduced in step (3), a single feature of 
a single individual is randomly chosen, and, with probability r, their corre- 
sponding trait is changed to a randomly selected value between 1 and q. 

In the absence of fluctuations, the system evolves towards absorbing 
states, i.e., frozen configurations that are not capable of further changes. 
However, for r > the system evolves continuously and, after a transient 
period, it attains a stationary state. Following previous studies on Axel- 
rod's model [2UI 125] , in this work we chiefly focus on systems of fixed size 
(N = 2500 nodes), fixed number of cultural attributes (F = 10) and fixed 
number of different cultural traits per attribute (q = 40). Furthermore, we 
also briefly discuss the effects (or lack thereof) observed by changing these 
model parameters. The results presented in the next Section correspond to 
observables measured over statistically-averaged ensembles in the station- 
ary regime, which were obtained by averaging over 200 different (randomly 
generated) initial configurations and 100 different network realizations. 

3 Results and discussion 

In order to set the stage for the investigation of modularity effects, let us 
first briefly summarize the main results concerning Axelrod's model defined 
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on the square lattice. As mentioned above, in the absence of fluctuations the 
system reaches absorbing configurations, in which the state of each individ- 
ual is fixed and not capable of further changes. With the inclusion of noise 
to model the effect of cultural drift, however, disordered multicultural con- 
figurations become metastable states that can eventually decay to ordered 
stable configurations [2S]. Whether this decay actually takes place or not 
depends on the competition between the noise rate, r, and the characteristic 
time for the relaxation of perturbations, T. For r < T _1 , repeated cycles of 
perturbation-relaxation processes drive the disordered system towards mono- 
cultural states, while, for r > T _1 , noise rates are large enough to hinder the 
relaxation mechanism, thus conserving the disorder. With arguments based 
on a mean-field description of a damage spreading process, the characteristic 
time for the relaxation of perturbations is estimated as T ~ iV IniV, where 
N is the system size [25] . 

In the absence of noise, the number of cultural traits is observed to play 
a key role in determining the final absorbing state: ordered monocultural 
configurations (for q < q c ) and disordered multicultural ones (for q > q c ) are 
separated by a finite critical value q c > [23, 25J. For instance, in a system 
of size N = 2500 with F — 10 cultural attributes, the transition takes place 
at q c ~ 50. For r > 0, however, the order-disorder transition solely depends 
on the effective noise rate r e ff = r x (1 — l/q). The very mild dependence 
on the parameter q just stems from the fact that, according to the third rule 
of the model dynamics, a perturbation can leave the cultural configuration 
unchanged with probability l/q. For the N = 2500 and F = 10 case, the 
order- disorder transition is observed around r e // ~ 10~ 4 [25J. 

When both noise and the mass media external field, M, are taken into 
account, interplay and competition effects are observed [22] • As the field in- 
tensity is increased, the transition shifts to higher noise levels. Since the or- 
dering is driven by the mass media field, the system attains a unique ordered 
state, namely, the monocultural state in which all individuals share their 
cultural traits with those of /i, the external message. In the absence of exter- 
nal fields, the noise-induced ordering leads to q F equally likely monocultural 
ground configurations, as well as to excursions from a ground configuration 
to another one. Note hence that, due to considerations of ergodicity and the 
multiplicity of ground states, it is not possible to simply define an "effective 
noise intensity" r' = r'(r, M) in order to trivially map the model with field 
onto an effective model without field. 

Let us now focus on the effects of modularity, which is modeled using a 
substrate topology that consists of two coupled random networks (see Sect. 2 
for details). Firstly, we will address cultural drift effects alone (r > 0, M = 0); 
later, we will study the case in which inhomogeneous mass media affect one 
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Figure 1: Typical snapshot configurations for different values of noise, r, 
and intercommunity connectedness, v, in the absence of mass media fields 
(M = 0). Nodes belonging to the community A (B) are shown on the left 
(right) side of each CRN realization. The most popular cultural state is 
shown in blue, the second most popular in green, the rest in greyscale. The 
network visualizations were created with Cytoscape [40J. 

of the communities (r > 0, M > 0) and explore the conditions for successful 
message spreading across the whole system. 

Figure 1 presents some typical snapshot configurations of the stationary 
regime for different values of noise, r, and intercommunity connectedness, 
v 1 in the absence of mass media fields: (a) r = 10~ 3 , u = 6 x 10 -2 ; (b) r = 
1CT 2 ,!/ = 4xl0~ 2 ; and (c) r = 10~ 3 , v = 1.5xlCT 2 . Here and throughout, the 
community A (B) is shown on the left (right) side of each network realization. 
For the sake of clarity, snapshot visualizations correspond to networks of 
small size (N = 100). The cultural state of the individuals is indicated 
by different node colors: the cultural state shared by the largest number of 
individuals is shown in blue, the second most popular cultural state is shown 
in green, while less frequent states are indicated in greyscale. 

The characteristic configuration for small noise levels and many intercom- 
munity links (Figure 1(a)) is a nearly full consensus: most of the individuals 
share the same cultural state. However, when considering larger noise rates 
(Figure 1(b)), the system undergoes a transition towards complete disorder, 
where the size of the most popular cultural state represents just a small frac- 
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tion of the total system's size. Indeed, these phenomena are reminiscent of 
the observed behavior of Axelrod's model in the lattice, where a finite crit- 
ical noise, r c , was found to separate the ordered, monocultural phase (for 
r < r c ) from the disordered, multicultural one (for r > r c ) [25]. However, 
strong effects arising from the modular structure of the substrate are ob- 
served at small values of u, leading to the appearance of a new phase (Figure 
1(c)). This is the ordered, bicultural phase, in which different cultural states 
prevail within each community. Interestingly, this behavior is in qualita- 
tive agreement with related work on a 2-state, majority rule model defined 
on substrates with community structure, where the coexistence of opposite 
opinions was observed [391 HI] • 

In order to quantitatively characterize different phases in the station- 
ary regime, we define A max (B max ) as the maximal number of members of 
community A (B) that share the same cultural state, normalized to unity. 
Furthermore, we define the vector a max (b max ) as the prevailing cultural 
state within community A (B). Once a stationary configuration is gen- 
erated, we classify the cultural state of each community as being ordered 
or disordered according to a simple majority criterion: A (B) is ordered if 
A m ax > 0.5 [B max > 0.5), and disordered otherwise. Moreover, when both 
communities are ordered, the whole system is in an ordered symmetric state 
if a max = b max , while it is in an ordered asymmetric state otherwise. Since A 
and B are indistinguishable communities, the combination of these different 
states leads to 4 possible phases. 

Figure 2 shows the resulting phase diagram in the r — v parameter space, 
in which the dominant phases for each region are displayed. At any given 
point on the r — v plane, the dominant phase is defined as the phase with 
the largest probability of occurrence (for instance, phase probability profiles 
for the M > case are shown in Figures 4 and 7 below). Thus, boundaries 
separating two phases correspond to states for which two dominant phases 
are equally probable, while triple points are associated to states for which 
three dominant phases are equally probable. The same definition was also 
adopted to determine the phase diagrams presented below (see Figures 5 and 
8). 

As anticipated, three distinct regimes prevail: a multicultural (disor- 
dered) phase, a monocultural phase in which both communities share the 
same cultural state, and a bicultural phase in which each community is or- 
dered, but in different states independent of each other. The effect of increas- 
ing r at a fixed value of v is that of increasing the disorder in the system, as 
expected. The noise-induced order-disorder transition is only mildly depen- 
dent on the number of links connecting both communities, with the transition 
curve located at r ~ 2 x 10~ 4 for v > 1.4 x 10~ 3 . The mild //-dependence of 
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Figure 2: Phase diagram for M = in the r — v parameter space. 

the order-disorder transition curve is due to finite-size effects: v plays here 
the role of tuning the "effective system size" from N e ff = N/2 (in the v —> 
limit, when the two communities are independent of each other) to N e ff = N 
(in the v —> 1 limit, when the modular structure washes out). 

It is within the region with small r where community structure effects 
become more noticeable. For small u, the two communities are weakly con- 
nected and do not influence each other, leading to the prevalence of the 
(ordered asymmetric) bicultural phase. However, increasing v the noise- 
driven ordered phases tend towards consensus, thus leading to a dominant 
(ordered symmetric) monocultural phase. A tri-critical point is found at 
(r = 2 x 10 -4 , v = 1.65 x 10 -3 ). As commented above, a qualitatively similar 
phase diagram was obtained in previous investigations of a 2-state majority 
rule model defined on substrates with community structure [39J, E] . 

The transition from the bicultural phase to the monocultural phase that 
takes place in the small-r region can be roughly estimated by the follow- 
ing theoretical argument. Under conditions of small rate of perturbations, 
we can assume that typically most of the nodes in the community A will 
tend to agree in the same cultural state a a (randomly chosen among any 
of the q F possible cultural vectors), and similarly, most of the nodes in the 
community B will share the cultural state <jb- The monocultural phase, 
°a = is driven by interactions between pairs of border nodes, i.e. those 
with intercommunity links. If a border node is chosen by rule (1) of the 
model dynamics, its interacting neighbor, chosen in turn by rule (2), has a 
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probability P ext = k ext /(k int + k ext ) to belong to a different community. If 
L u is the total number of intercommunity links, we can assume that the in- 
teraction between communities is effectively present for P ext L v > 1. Thus, 
PextLy ~ 1 can be taken as a rough estimate for the occurrence of the mono- 
cultural/bicultural phase transition in the small-r region. Since, on average, 
kext = vk int and v <C 1, a border node has k ext = 1, i.e. only one external 
neighbor. Using L u = 2uN, the condition P ext L u ~ 1 reads 



Replacing k int = 4 and N = 2500, we obtain v ~ 10~ 3 , which provides a good 
estimate for the boundary observed in Figure 2 between the monocultural 
and bicultural phases. An immediate consequence of this simple theoretical 
argument is that results for different network sizes should scale with v and 
N through L v oc vN . This predicted behavior was indeed confirmed by our 
simulations of systems of different size. Moreover, increasing N also shifts 
the boundary between ordered phases and the multicultural phase towards 
smaller values of noise, which is consistent with previous observations of 
noise-driven transitions in Axelrod's model defined on the square lattice [25] . 

Let us now address the case in which mass media affect one of the commu- 
nities and explore the conditions for successful message spreading across the 
whole system. In order to capture the characteristic behavior of this system 
in the multidimensional parameter space, we consider separately the small- 
r ordered region and the large-r disordered case. Since r can be regarded 
as a measure of the intrinsic individual determination or "free will" relative 
to the influence exerted by neighbors and mass media, the small-r scenario 
represents a society with individuals subject to strong social pressure, while 
the large-r case corresponds to a society characterized by loose social ties. 
Within these different scenarios, the adoption of inhomogeneous mass media 
fields that introduce a physical distinction between the dynamics of both 
communities will drive the system across different phase transitions. 

Figure 3 shows typical snapshot configurations of size N = 100 for dif- 
ferent values of M and v within the small-r regime (r = 10~ 3 ): (a) M = 
KT 3 , v = 6 x 10~ 2 ; (b) M = 10~ 2 , v = 6 x HT 2 ; (c) M = 10" 3 , u = 2 x KT 2 ; 
and (d) M = 10~ 2 , u = 2 x 10 -2 . Recall that nodes belonging to the commu- 
nity A (B) are shown on the left (right) side of each network realization. The 
cultural state that corresponds to the external message, /i, is shown in red. 
Other states are shown in blue (most popular among non-/z states), green 
(second most popular) and greyscale (all other states). 

When communities are strongly interconnected, the system evolves to- 
wards consensus, where most of the individuals share the same cultural state 
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Figure 3: Typical snapshot configurations for different values of M and u 
within the small-r regime. The //-state is shown in red, the most popular 
among non-/x states appears in blue, the second most popular non-/x state in 
green, while other states are shown in greyscale. 



(Figures 3(a)-(b)). However, the nature of the attained consensus depends 
on the strength of the mass media field: for small M, the system organizes 
itself into any of the q F possible monocultural states, while increasing M a 
transition takes place towards a regime dominated by the mass media mes- 
sage. When communities are instead sparsely interconnected, they tend to 
evolve independently of each other (Figures 3(c)-(d)). However, analogously 
to the case where communities are tightly bound together, the system under- 
goes an M-driven transition from a regime where communities are ordered 
but independently organized into different non-/x states (Figure 3(c)) to a 
phase where the mass media message prevails within community A, while 
the community B is in a non-/x state (Figure 3(d)). Notice that, lacking 
enough intercommunity links, even a very intense mass media campaign will 
fail to convey its message to the whole society. 

In order to quantitatively distinguish among different phases, we can 
follow a procedure similar to that described above for the M = case. 
However, additional phases now arise from the fact that the cultural state 
corresponding to the mass media message, /x, is physically distinguishable 
from the other q F — 1 possible cultural states. For instance, the community 
A can be either dominated by the mass media message (A max > 0.5 and 
a max = fJ>), ordered in a different cultural state (A max > 0.5 and a max ^ 
fx), or disordered (A max < 0.5). Taking also into account the distinction 
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Figure 4: Probability of occurrence of the relevant phases as a function of M 
for the small-r regime (with r = 10~ 5 ) and different values of the community 
connectivity: (a) v — 6 x 10~ 4 and (b) v — 2 x 10~ 3 . 



between symmetric and asymmetric ordered states (which is relevant when 
A and B are ordered in cultural states both different from /z), this ultimately 
leads to 10 possible different phases. However, as suggested by the snapshot 
configurations shown in Figure 3, only 4 phases are relevant. 

Figure 4 shows the probability of occurrence of the relevant phases as 
a function of the message intensity, corresponding to the small-r regime 
(with r = 10~ 5 ) and for different values of the community connectivity: 
(a) v — 6 x 10~ 4 and (b) v — 2 x 10~ 3 . In agreement with our qualitative 
discussion, Figure 4(a) shows that different kinds of asymmetric phase pre- 
vail when communities are loosely interconnected. Indeed, for small mass 
media fields, communities are predominantly in different non-/i states, i.e. 
the (A„ on _ M , -B„ „- m )a phase, while increasing M above M c = 7 x 10~ 4 the 
system is most often found in a bicultural phase with fj, prevailing within 
community A, labeled as the (A^, B non ^^) phase. Note also that, somewhat 
counterintuitively, the probability of achieving overall consensus tends to de- 
crease as a function of the mass media intensity. This phenomenon is due to 
the fact that the external field prevents the independent auto-organization of 
the whole system in a non-/i state, while the lack of strong intercommunity 
ties prevent the message from reaching out far beyond the region directly 
exposed to the inhomogeneous mass media field. Figure 4(b) shows that, on 
the contrary, strongly interconnected communities allow the whole system 
to reach consensus. Increasing the mass media field, the system undergoes 
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Figure 5: M — v phase diagram within the small-r regime (for r = 10 5 ). 

the expected transition from non-/x monocultural states, i.e. the symmetric 
(v4 no „_ M , B non -^) s phase, to /^-consensus, indicated as (A^B^). 

The corresponding phase diagram in the M — v parameter space is shown 
in Figure 5. In the low-M end, we observe that, increasing the intercommu- 
nity connectedness, the dominating phase changes from asymmetric non-/x 
to symmetric non-yU. In fact, this transition matches the bicultural to mono- 
cultural phase transition observed earlier in the small-r end in the absence of 
mass media fields (recall Figure 2). As discussed above, successful message 
spreading across the whole system can only be achieved when sufficiently 
strong mass media fields are applied on a sufficiently interconnected system. 
The boundaries for the //-consensus region are approximately M > 10~ 3 and 
v > 1.5 x 10~ 3 . These results stress the fact that, in a general scenario, 
well-designed, cost-effective advertising campaigns should take into account 
the specific modular structure of the target population. Indeed, even very in- 
tense (and, hence, costly) mass media campaigns may fail if social modularity 
effects are disregarded. 

Considering networks of different size, phase diagrams can be made roughly 
invariant along the vertical axis by adopting the scaling relation uN . Indeed, 
this indicates that the relevant quantity defining the actual degree of inter- 
connectedness is the total number of intercommunity links, L„, in agreement 
with the theoretical argument presented above, Eq.flT]). The boundary be- 
tween non-/i states and the \x— consensus phase, which in Figure 5 appears 
around M ~ 1CT 3 , is observed to shift towards lower values of M as the 
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among non-/x states appears in blue, the second most popular non-/z state in 
green, while other states are shown in greyscale. 



network size is increased. We also explored the stability of our results under 
changes in the number of cultural attributes, F, and the number of different 
cultural traits per attribute, q, without noticing any significant variations. 
This behavior agrees well with previous investigations on Axelrod's model, 
where results roughly independent of the parameter F [23] and the parame- 
ter q (for the model with noise, provided that g > 1) [251 [29] were reported. 
A similar behavior is also observed in the large-r regime, which is discussed 
below. 

As anticipated, we now consider the large-r scenario, which corresponds 
to a society characterized by loose social ties. Figure 6 shows typical snapshot 
configurations for different values of M and v within the large-r regime (r = 
1(T 2 ), with the same coloring scheme used in the visualizations of Figure 3. 
The corresponding parameter values are: (a) M = 1.5 x 10~ 2 , v = 6 x 10~ 2 ; 
(b) M = 10-\u = 6 x 10" 2 ; (c) M = 1.5 x 10~ 2 , is = 2 x 10~ 2 ; and (d) 
M = lO-\v = 2 x 1(T 2 . 

Irrespective of connectivity, the low-M region is characterized by disor- 
dered multicultural configurations (Figures 6(a) and 6(c)). Indeed, disor- 
dered states are characteristic of the large-r region in the absence of mass 
media fields (compare to Figure 1(b)). Within this scenario, order can be 
achieved only when strong external fields oppose the large intrinsic noise (see 
Figure 6(b)). However, if the communities are sparsely interconnected, even 
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Figure 7: Probability of occurrence of the relevant phases as a function of M 
for the large-r regime (with r = 10" 3 ) and different values of the community 
connectivity: (a) v — 2 x 1(T 3 and (b) v — 1.7 x 10~ 2 . 



strong fields are not capable of driving the system towards consensus: while 
the mass media message prevails within community A, the community B is 
instead in a disordered multicultural state (Figure 6(d)). 

Following the procedure described above, we can determine the probabil- 
ity of occurrence of each phase as a function of the field intensity. However, 
out of 10 possible phases, only 3 are relevant. These are shown in Fig- 
ure 7 for r = 10 -3 and different values of the community connectivity: (a) 
v = 2 x 10~ 3 and (b) v — 1.7 x 10~ 2 . As discussed above, when the com- 
munities are loosely interconnected (Figure 7(a)) the system is intrinsically 
disordered. Only a strong field can oppose the large noise level and order the 
community A, thus leading to a transition from the (Adi S , B<u s ) phase to the 
{A^Bdis) phase. For highly connected communities (Figure 7(b)), instead, 
the strong field is able to order the whole system in the //-state, driving the 
phase transition from (A dis ,B dis ) to (A^B^). 

Figure 8 shows the phase diagram in the M — v parameter space cor- 
responding to the large-r regime (r = 1CT 3 ). Matching the large-r region 
for the M = phase diagram of Figure 2, the low-M end is dominated by 
the multicultural disordered phase. Increasing M, two different phases can 
be reached depending on the intercommunity connectedness: the (A^, Bdi S ) 
phase, for loosely connected communities, and the //-consensus, when the 
communities are more strongly bound together. 

Finally, let us discuss some subtle, intriguing effects that result from the 
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Figure 8: M — v phase diagram within the large-r regime (for r = 10 3 ). 



model dynamics. Recalling the phase probability distributions of Figure 4(b), 
the plot corresponding to the (A^, B^) phase shows a dip at large message in- 
tensities, which is correlated with the occurrence of a bump in the plot of the 
(A^, B non -n) phase. Far from being an artifactual feature due to poor statis- 
tics, this phenomenon stems from the dynamical rules of the model and can 
be understood on the basis of a sound sociological interpretation. According 
to the second dynamical rule, within the community A, the parameter M 
regulates the competition between two different types of interaction: that 
of an individual with the mass media, and that between neighboring indi- 
viduals. This implies that, besides the ordering effect driven by the mass 
media interaction, which tends to align all cultural traits with the external 
message, there is also a competing disordering mechanism: individuals sub- 
ject to strong mass media fields have a low probability of interacting (and, 
thus, of increasing the similarity) with their social neighbors. Although the 
former tends to prevail and is, ultimately, the mechanism responsible for the 
yU-consensus ordering observed in the large-M end of the phase diagrams, 
competition effects lead to visible features such as the dip and the bump 
noted above. 

In order to confirm this explanation, we computed the effective link 
weight, w e ff(A — B), as the mean cultural overlap between neighbors that 
belong to different communities, and compared it to the effective link weight 
within each community. Figure 9 shows the effective link weights as a func- 
tion of M for the small-r regime (with r = 1CT 5 ) and a large connectiv- 
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Figure 9: Effective link weights as a function of M for r = 10~ 5 and 
v = 2 x 10~ 3 (solid lines). Link weight plots are shown separately for in- 
tracommunity connections within each community, as well as for intercom- 
munity connections. For comparison, the probability of occurrence of the 
relevant phases are also shown (dashed lines). See more details in the text. 



ity {y — 2 x 10~ 3 ). For the sake of comparison, the probability of occur- 
rence of the relevant phases (same as those in Figure 4(b)) are also shown. 
Marked with arrows, we indicate two distinct features in the plots of effec- 
tive link weight: a first dip (local minimum) in all three plots taking place 
at M — 10~ 3 , and a second dip observed at M — 2.6 x 10~ 3 only in the 
plot of intercommunity links. The former is well correlated with the phase 
transition from {^A non _^ B n(m _^)s to (A^, B^), and hence reflects the corre- 
sponding phase changes within each community. Instead, the latter is well 
correlated with the dip in the probability of the (A^, B^) phase, as well as 
with the bump in the probability of the (A^, £> non _ M ) phase. 

These results show that individuals subject to intense mass media fields 
are less likely to interact with their social neighbors, hampering message 
spreading processes across community boundaries. Thus, induced by mass 
media pervasiveness, the tendency of individuals towards isolated behavior 
is captured and well accounted for by this model. 
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4 Conclusions 



In the context of an extension of Axelrod's model for social influence, we 
studied how the modular structure of social networks affects the propagation 
of mass media messages across the society. The community structure of social 
networks was represented by coupled random networks, in which two random 
graphs are connected by intercommunity links. 

In the absence of mass media, we observed the prevalence of three distinct 
phases, depending on the values of cultural drift (i.e. the level of intrinsic 
noise) and community interconnectedness: the ordered monocultural phase, 
the ordered bicultural phase, and the disordered multicultural phase. The 
obtained phase diagram is qualitatively similar to that reported for the ma- 
jority rule model defined on modular substrates (39j SI]. 

Then, considering inhomogeneous advertising campaigns, we studied the 
system's response to the spreading of the mass media message. We considered 
separately two different scenarios: the small noise regime, which represents 
a society with individuals subject to strong social pressure, and the large 
noise regime, which is characterized by loose social ties. Incorporating the 
intensity of the mass media field as an additional parameter, we observed the 
emergence of new phases, which led to a very rich, novel phase diagram in 
the multidimensional parameter space. Our results show that the design of 
successful, cost-effective advertising campaigns should take into account the 
specific modular structure of the target population. Indeed, even very intense 
(and, hence, costly) mass media campaigns may fail if social modularity 
effects are disregarded. 

Certainly, the simplified scenario considered here leaves room for fur- 
ther investigations that may look for modularity effects when communities 
of different size and different inter- /intra-community connectedness are con- 
sidered. In this vein, the study of social influence and message spreading 
on real modular substrates taken from social interaction data (such as e.g. 
large-scale mobile phone usage with detailed space-time resolution [42, 43j) 
could lead to interesting new results. We thus hope that the present findings 
will contribute to the growing interdisciplinary efforts in the mathematical 
modeling of social dynamics phenomena, and stimulate further work. 
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